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Q ' We present a detailed calculation of channeling radiation of planar-channeled 

positrons from crystal targets in the framework of our approach, which was pro- 
posed recently. In contrast to previous calculations of channeling radiation in crys- 
tals, our calculation takes into account the interference between different transition 
amplitudes. The development stemmed from the idea that the amplitude for a given 
process is the sum of the transition amplitudes for each transition to lower state 
I of transverse energy with the same energy differences between bound-bound tran- 

' sitions. It seems that a consistent interpretation is only possible if positrons move 

in a nearly harmonic planar potential with equidistant energy levels. 
I Key words: Quantum theory of channeling; Radiation from positrons; Har- 

monic potential 
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=5 . 1 Introduction 

Charged particles directed into a crystal approximately parallel to one of the crystal 
planes will be planar channeled. For positively charged particles, such as positrons, the 
r> I channel is between the crystal planes. In classical terms, the particle's momentum forms 
I some small angle 6 with respect to the crystal plane. This angle must be less than the 
Lindhard's critical angle 6p for planar channeling to occur. The motion of the particles then 
consists of a periodic back-and-forth reflection of the boundaries of the planar channel. 
An accelerated periodic motion of this kind will lead to the emission of radiation. From a 
quantum- mechanical viewpoint, the channel is the source of a one- dimensional potential 
well for planar channeling in the direction transverse to the particle's motion, which 
gives rise to transversely bound states for the particle. Its transitions to lower levels are 
accompanied by the emission of radiation with frequencies related to the energy differences 
of the levels. Theoretical studies of the radiation from planar-channeled electrons and 
positrons are due to M.A. Kumakhov [1,2], N.K. Zhevago [3], A.I. Akhiezer, LA. Akhiezer 
and N.F. Shul'ga [4,6], A.W. Saenz, H. Uberall and A. Nagl [5], and to V.A. Bazylev, 
V.V. Beloshitsky, V.I. Glebov, N.K. Zhevago, M.A. Kumakhov and Ch. Trikalinos [7]. 
Experimentally, the channeling radiation of positrons was observed by the different groups 
[8-11], demonstrating strong and sharp peaks in the spectrum. The purpose of the present 
work is to calculate the spectral-angular distribution of the channeling radiation intensity 
emitted from positrons in the framework of approach which was proposed recently [12]. 

* Corresponding author. E-mail address: mshatnev@yahoo.com. 
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2 Wave functions and transverse potential 



The following quantum mechanical calculation of planar positron channeling radiation 
utilizes the developments of Kumakhov and Wedell [2] , and Zhevago [3] . One begins with 
the time-independent Dirac equation for a relativistic particle moving with momentum 
p II = {0,Py,Pz) in a one-dimensional planar potential V(x) periodic in the x-direction 
(which is normal to the channeling planes) 

• V + - /3m)^' = (1) 

where m and E arc the particle's mass and energy, 'a and (3 are the Dirac matrices. 
Separating the wave function into large and small components, 



and using the standard representation for the Dirac matrices, leads to a Pauli-type equa- 
tion for the large components, 

'a ■ V{E - V{x) + m)-^^ ■ V^a + {E - V{x) - m)^^ = 0. (3) 

Since a potential V{x) is independent of y and z, the solution of Eq. (3) is a plane wave 
in the yz plane 

oc exp{i{pzZ + pyy))(p{x)x. (4) 

We now take advantage of the fact that the particle energy is £^ > 1 GeV, which is much 
larger than the planar potential energy, which is on the order of 10 eV. This allows us 
to transform Eq. (3) into a one-dimensional, relativistic Schrodinger equation for the 
transverse motion 

"2^^^ + ^^'^^'^^'^^^"'^^'^^' 

where 

E"^ -m^ - pi- pI , ^ 



In a given potential, the latter will assume certain bound-state eigenvalues 6„ < (n = 
0, 1, 2...), with corresponding eigenfunctions (^nl^:^)- While Eq. (6) then seems also to lead 
to a quantization of = ^p\ + p^, this is not the case in the approximation described 
after Eq. (4) where e <^ E; here, Eq.(6) shows that to sufficient accuracy. 



E^E^^ = ^J^\^ + w?. (7) 
The wave function of Eq. (2) is finally obtained in the form 

J exp(ip||-r||)9?„(a;), (8) 



where 



being the two-dimensional normalization volume for the plane waves of Eq.(4), and x 
being a two-component spinor which is (q) or (°) when the particle spin points in the +z 
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or in the —z direction in the rest frame, respectively. Since V(x) is a periodic potential, 
we may consider the problem in a single well with harmonic potential Vix) = Vq-x^ which 
describes planar channeling of positrons along the (110) plane of a silicon crystal [2,3]. 
For this potential, it is well known that Eq. (5) can be solved, and the corresponding 
eigenfunctions being given by 

<Pn{x) = \ —^^e^v{-E^x'|2)Hr^{^^mx), (10) 
V TT V2"n! 

where are the Hermite polynomials, Q is the oscillation frequency given by 

where dp is the distance between planes in the corresponding units, and Vq is the depth 
of the potential well. The corresponding eigenvalues being given by 

£„ = fi(n + l/2). (12) 



3 Matrix element and intensity of channeling radia- 
tion 

The crystal plane which channels the positrons is the yz plane, and the longitudinal 
component p\\ of the positron momentum is taken parallel to the ;i;-axis {py = 0). The 

emitted photon k makes an angle 9 with p^\, and the azimuth between the emission 

( ^ ■ plane and the crystal plane is (f. Following Kumakhov [2], we can choose two 
linear polarization vectors of the photon in the next way:T^i being normal to the emission 
plane and ~e2 lying in the emission plane. The differential emitted photon intensity for a 
given polarization is, by Fermi's golden rule, 

V(Pk L'^d'^p'n 

d'h = 27ru;\M,f\'5{u;-u;(3\\COse- u;^^,)^"^, (13) 

where V = is the normalization volume, A = 1, 2 indicates the linear photon polariza- 
tion, uj is the frequency of the radiation, k is the photon momentum, ^„„'= Sn — £n' 
the transition energy between the two levels n and n in the rest frame of the positron, 
and p II = {0,Py,p^) are the components of the final positron momentum vector parallel 
to the crystal planes. The transition matrix element is given by [13] 

|M.r^^|//. (14) 

where — and 

Jx^ j ^'^ctAe-^'^'^^dV, (15) 
with ax = 'a ■ ~s'*y The prime indicates the final state. Momentum conservation 

k = P^ii - 'p'w, (16) 



3 



and energy conservation, 

^1, +£n = ^li+£n'+^, (17) 

lead to the "Doppler shift" formula for the photon energy uj^^' corresponding to the 
transition n ^ n , 



nn 



(18) 



l-/?||COS^' 

II 

where j3\\ = in the approximation o;^^' <S E. This approximation is justified for the 

II 

case E ~ GeV and £n — < 30 considered here. Since fti (1 — 

(where 7 = :^), cos ^ (1 — |^^) and — = Q{n — n) Eq. (18) can be expressed as 

It follows from Eq. (19) that the radiation of a maximum frequency 

a;„„,=272Q(n-n') (20) 

is emitted in the forward direction (at = 0). The case n — n = 1 corresponds to the peak 
values of the experimental channeling radiation spectra [11], being the first harmonic with 
the photon energy u; — As it follows from Eq.(19), photons emitted via positron 

transition from any initial level n to the final level n — 1 are identical (i.e., have the same 
energies for the same emission angles). This means that the resulting amplitude should 
be given by an additive superposition of amplitudes of all such transitions. The positron 
state outside the crystal {z < 0) is a plane wave, whereas inside the crystal {z > 0), the 
part of its wave function corresponding to the transverse motion is a superposition of the 
harmonic oscillator eigenvectors. Factors c„ describing transitions from the initial state 
to states with the transverse energy levels n can be found using boundary conditions set 
upon the wave function at the crystal boundary {z = 0). Then, a transition to the closest 
lower level n — 1 occurs with emission of a photon having energy uo. One may expect that 
the total amplitude of the transition from the initial to final state accompanied by the 
photon emission is determined by products of the amplitudes c„ and M„_„_i. Following 
the rules of the quantum mechanics, we express this amplitude as 

A(xJ2cnMn,n-l: (21) 
n 

were summation is performed over all the harmonic oscillator levels. Also we must find 
an additive superposition of amplitudes of all transitions n — > n — 2, n — 3, ...etc. Taking 
into account these considerations, we can write the transition matrix element in the form 



^ ^ Cn-^n, n—'j 



(22) 



where j — n — n . According to Eq. (15), we find, using Eq. (8) for the wave functions, 
the first-order matrix element corresponding to the n — > n — j transverse transition 

Ja = (27r)25(^|| - tn - h)NN'x*l'l(A + i[B^])x, (23) 
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where 



Bx 



•7-(2) 



1 



1 



AD ^ , 1 1 
1 ■ V 

n,n—j r 



E + m E — uj + m ^ 
) + 



r(i) 



n,n-j \ 



P 



P 



■' E -\-m E — uj + m 



with 



n,n—j 

r(2) 
n,n-j 



J exp{-ikx:x)ipl_j{x)(pn{x)dx, 
exp{-ikxx)(p*^_j (x) ^'^J^^^ dx, 



(24) 
(25) 
(26) 

(27) 

(28) 
(29) 



Then we calculate the matrix element and the differential intensity, and after summing 
over polarization of emitted photon and the positron using the formulae 



X X 

X X 



Sik - riiUk, 
1 + ^^ ^ 1 
2 ^2' 
1 + ^T' 



(30) 
(31) 

(32) 



where ~^ and arc the unit spin vectors of the initial and final positron (in the rest 
frame of the positron), respectively, and taking into account that '^nn'~ 3^ according to 
Eq. (12), we find by integrating over d?p'» 



(PI 



dudo 27r 



E 



d{u! — a;/3|| cos 6' — jfl), 



(33) 



where the factors c„ are, in the case of the parabolic potential and when the initial positron 
is a plane wave, given by 



V2"-in! V Efl 



(34) 



4 Conclusion 

Following Kumakhov and Wedell [2], and Zhevago [3], the spectral-angular distribution 
of emitted photons is represented as 



d^I 
dujdo 



oc 



(35) 



The sum entering Eq. (35) is the one over the quantum numbers / of the transverse 
motion of the particle. Then, the probability of having a definite transverse energy is 
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taken into account by multiplying each term of this sum by a corresponding factor. In 
our consideration, discrete levels of the transverse motion in the harmonic oscillator po- 
tential refer to the intermediate state of the particle. Accordingly, the contribution to the 
intensity due to transitions, e.g., between the closest levels is determined by the square 
of the absolute value of Eq. (21) 



dujdo 



(36) 



In other words, unlike Ref. [2,3], we get an expression that contains interference terms 
mixing amplitudes of photon emission from different equidistant levels. We would like to 
note that dynamics of channeling electron in a crystal differs from that of the positron 
case. The transverse potential well for the electron does not give rise to equidistant energy 
levels for transverse particle motion. Therefore, there are no interference contributions 
to the photon emission intensity similar to those present in Eq. (36). In our opinion, 
this could explain the greater intensity in case of channeling positron compared to that 
for the electron observed in experiment [11]. We think that studying the problems like 
the one under consideration in this paper can help in the development of methods for 
obtaining polarized photon beams. Corresponding numerical calculations will be given in 
a subsequent publication. 
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